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Abstract. Extensive simulations are performed to study the persistence behavior of a conserved lattice gas 
model exhibiting an absorbing phase transition from an active phase into an inactive phase. Both the global 
and the local persistence exponents are determined in two and higher dimensions. The local persistence 
exponent obeys a scaling relation involving the order parameter exponent of the absorbing phase transition. 
Furthermore we observe that the global persistence exponent exceeds its local counterpart in all dimensions 
in contrast to the known persistence behavior in reversible phase transitions. 
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1 Introduction 

In the last decade numerous articles have been devoted 
to the study of persistence behavior in nonequilibrium 
systems |],||,|,|p,|,|,[o[|l3. The basic quantity of 
interest is the persistence distribution P{t), defined as 
the probability that a certain physical quantity does not 
change its state up to time t during a stochastic evolu- 
tion. The persistence distribution has been observed to 
decay as P{t) ^ t~^, where 9 is denoted as the persis- 
tence exponent. Since the persistence behavior depends 
on the entire history of evolution of the nonequilibrium 
process, there are only few instances of an exact calcula- 
tion of the persistence exponent in the hterature In 
general the persistence exponent is nontrivial in the sense 
that it cannot be obtained from scaling relations, except of 
certain Markovian systems jij. There exists two different 
types of persistence exponent 0i and 6^ associated with 
the persistence behavior of the local and the global order 
parameter, respectively. These two exponents have been 
observed in general to be different from each other. Usu- 
ally for reversible phase transitions the inequality 9g < 9\ 
is observed, however in the case of the continuous absorb- 
ing phase transition of the Ziff-Gulari-Barshad model jl^ , 
which belongs to the universality class of directed perco- 
lation^ the reverse inequality has been observed numeri- 
cally 0. 

Recently Rossi et al. introduced a conserved lattice 
gas (CLG) model that exhibits a continuous phase transi- 
tion to an absorbing state at a critical value of the parti- 
cle density |l3t . Numerical simulations in various dimen- 
sions 10 confirm the conjecture of |l3| that the CLG 
model belongs to a new universality class different from 



the well known universality class of directed percolation 
(see for an overview). In this work we consider the 
persistence distributions of the CLG model in various di- 
mensions {D = 2,3,4,5). We determine numerically the 
global and local persistence exponents as well as the ex- 
ponent I'll which describes the divergence of the correla- 
tion time at the transition point. To our knowledge this 
is the first time that the persistence behavior of a non- 
equilibrium critical system with absorbing states is inves- 
tigated in higher dimensions, i.e., below and above the 
upper critical dimension D^- Similar to we observe 
the reverse inequahty 9g > 9\ for the CLG model. We be- 
lieve that systems exhibiting irreversible phase transitions 
into an absorbing state are generally characterized by this 
reversed inequality in contrast to systems exhibiting a re- 
versible phase transition, in particular equilibrium phase 
transitions. 

The paper is organized as follows. In the next section 
we describe the CLG model. In section ||we define the lo- 
cal and the global persistence distributions for this model 
and describe their steady state scaling behavior in the crit- 
ical regime. In the following two sections we present the 
determination of the global (sec. ^) and local (sec. ||) per- 
sistence exponents. With a short discussion we conclude 
the paper in section ^. 



2 The CLG model 

We consider the CLG model (see [^) on D-dimensional 
cubic lattices of linear size L. Initially one distributes ran- 
domly N = pL particles on the system where p denotes 
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the particle density. In order to mimic a repulsive interac- 
tion a given particle is considered as active if at least one 
of its 2D neighboring sites on the cubic lattice is occupied 
by another particle. If all neighboring sites are empty the 
particle remains inactive. Active particles are moved in the 
next update step to one of their empty nearest neighbor 
sites, selected at random. Starting from a random distri- 
bution of particles the system reaches after a transient 
regime a steady state which is characterized by the den- 
sity of active sites p^. The average density (pa) is the order 
parameter of the absorbing phase transition, i.e., it van- 
ishes if the control parameter p is lower than the critical 
value pc- Approaching the critical value from above the 
order parameter scales as (pa) ~ Sp^ , where 6p — p — pc- 
The order parameter exponent (3 as well as the scaling be- 
havior of the order parameter fluctuations are investigated 
numerically in various dimensions in [Q. In our simula- 
tions, a random sequential update scheme is applied, i.e., 
all active sites are updated at each time step in a ran- 
domly chosen sequence. The considered system sizes are 
L < 1024 for D = 2, L < 128 {D = 3), L < 32 {D = 4) 
and L < 16 in the case oi D = 5. All measurements are 
performed in the steady state and the used equilibration 
procedures are discussed in detail in . 

3 Persistence distributions 

In this work we consider the persistence behavior of the 
CLG model. Usually one distinguishes the so-called global 
persistence distribution Pg{t) and the local persistence 
distributions Pi{t). The global persistence characterizes 
the time evolution of the order parameter, i.e., it char- 
acterizes how long the whole system remains in a given 
phase. For instance in magnetic systems Pg(t) can be de- 
fined as the probability that the magnetization changes 
sign after the duration t. At criticahty Pg{t) is expected 
(see for instance Q) to decay algebraically 

p,{t) - t-'s (1) 

where the time t corresponds to the number of update 
steps. Off criticahty, the finite correlation time limits 
this power-law behavior and the corresponding probability 
distribution displays a cut-off at t = C(Cii). Using ~ 
(5p~''ii and standard scaling arguments one expects that 
the persistence probability obeys the ansatz 

Pg(t) = ^p''"^- PgiSp"^^ t), (2) 

with the universal function Vg. Thus the scaling behavior 
of the probability distribution Pg{t) allows to determine 
beside the persistence exponent 9g the exponent of the 
temporal correlations in all considered dimensions. 

The phase transition of the CLG model is an irre- 
versible transition into an absorbing state, i.e., the or- 
der parameter (pa) is a strictly positive quantity. In this 
cases it is usual [|o[|l^ to measure the probability that 
the deviation of the density of active sites from the order 
parameter, i.e. pa— (pa), changes sign after t update steps. 



Furthermore, the distribution of the deviations could be 
different below and above (pa). In order to take this asym- 
metry into account we calculate two different quantities 
Pg{t) and Pg{t) corresponding to the global persistence 
distribution for Pa(i) > (pa) and Pa(t) < (pa), respectively. 

The so-called local persistence is defined as the prob- 
ability that the local order parameter at a given site, for 
instance a single spin, remains in a phase, i.e., it does not 
flip within a duration t. In the case of the CLG model 
we consider the steady state probability distribution Pi(t) 
that a particle remains inactive for exactly t update steps. 
Since no spontaneous generation of active sites is allowed 
in the CLG model, a given site may remain inactive for a 
long time until it is re-activated in the presence of another 
active site. Thus one expects that the corresponding prob- 
ability distribution decays algebraically with an exponent 

Pi{t) t-'^. (3) 

Alternatively one may ask how long a given particle re- 
mains active. But active sites can be spontaneously "anni- 
hilated" to become inactive and the corresponding prob- 
abilit y d istribution decays exponentially. This "asymme- 
try" in the dynamical behavior of active and inactive 
sites is already known from other systems exhibiting an 
absorbing phase transition, for instance directed percola- 
tion. 

Above the transition point (p > pc) the local persis- 
tence distribution displays a cut-off at t — 0{£,\), where 
^1 denotes the correlation time of the effective single site 
process. Assuming that scales as ^ 6p~'^' one obtains 
the scaling behavior of the local persistence distribution 

Pi{t) - Sp-^^'^ ViiSp""' t) (4) 

with the local critical exponents 6i and vi. These local ex- 
ponents are connected to the order parameter exponent p. 
Using the scaling ansatz Eq. one can calculate the av- 
erage time between two activation events 

(t) = J tPi{t)dt ~ Jp-'^'^^-oo. (5) 

Approaching the transition point {Sp — > 0) the average 
time between two activation events diverges (we will see 
that 01 < 2 in all dimensions). The number of activation 
events within a certain duration is proportional to 1/ (t) . 
On the other hand the number of activations is propor- 
tional to the probability that a site is active at a given 
time. Therefore l/{t) oc (pa) and this leads to the scaling 
relation 

i/i (2 - 9i) = (3. (6) 

In spatially extended systems the time evolution of a 
single site is coupled to that of its neighbors, thus the 
effective single site evolution is non-Markovian. In this 
sense the local persistence distribution is related to an in- 
finite point correlation function (see for instance pT|,p^) 
whereas the persistence distribution of the order param- 
eter is related to a two-point correlation function. Thus 
one expects that the exponent vi does not agree with the 
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exponent of the order parameter i^i ^ i/y . But this in- 
equahty does not imply that the order parameter scahng 
behavior of the CLG model is characterized by an addi- 
tional second correlation time because vi just describes the 
temporal correlations of an effective single site time-series 
and not of the temporal evolution of the order param- 
eter. Furthermore, one may expect for the same reason 
that the persistence exponents are different, i.e., 6g ^ 9i. 
Previous investigations suggest a general behavior with 
9g < 9\. In this case the persistence behavior of the global 
properties are stronger than those of the local proper- 
ties. For instance the average duration between two spin 
flips in ferromagnetic systems is smaller than the dura- 
tion between order parameter flips. Recently, a violation 
of this inequality was reported by Hinrichsen and Kodu- 
vely who observed 9g « 9\ in the case of directed perco- 
lation jl^. Additionally the persistence distributions of 
the Ziff-Gulari-Barshad model which exhibits a con- 
tinuous absorbing phase transition belonging to the di- 
rected percolation universality class, are characterized in 
the "long-time regime" by the inequality 9g > 9i for 
D = 2. 




t (P-Pc)'« 

Fig. 2. Scaling analysis of the two-dimensional global persis- 
tence distributions P^{t) and Pg{t). In the case of P^{t) the 
curves are shifted in the lower-left direction. The values of the 
exponents are presented in the text. 



4 Global persistence distributions 

We consider the global persistence distributions for dimen- 
sions D = 2, 3, 4. As explained above one expects a cut-off 
of the probability distributions for t w 0(^11)- Thus, ©(^y) 
update steps have to be performed to measure one of these 
long-time events. This shows that extensive simulations 
are required to obtain the global probability distribution 
with a sufficient accuracy. For instance 5 x 10^ single up- 
date steps were used in to determine the critical be- 
havior of the order parameter (pa) for D = 2. Here, we 
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Fig. 1. The global persistence distributions Pg{t) and Pg{t) 
for a given system size and particle density. The different asym- 
metric behavior below and above the average value (pa) (see 
text) can be clearly seen in the cut-off behavior. 



used for similar system parameters 4 x 10^ update steps 
in the steady state which yield approximately 6 x 10^ sin- 
gle events for the measurement of the global probability 
distribution. 

Fig. [l] shows the distributions Pg{t) and Pg{t) for the 
two-dimensional model in a log- log plot. As one can see the 
slopes of the curves are nearly the same but the cuf-off be- 
havior differs significantly. Simulations for different values 
of 5p are performed and the obtained curves are analyzed 
according to the scaling ansatz Eq. (j^). In the case of the 
two-dimensional model we obtain a good data collapse for 
91 = 1.57 ± 0.05, i^if = 1.15 ± 0.04 and 9< = 1.60 ± 0.06, 

= 1.18 ± 0.05, respectively. The corresponding curves 
are shown in Fig.^. Although the distributions Pgit) and 
Pg{t) differ at the cut-offs the scaling behavior, i.e. the 
exponents, of both distributions agree within the error- 
bars. 

Our simulations yield that the different cut-off behav- 
ior of Pg{t) and Pgit) vanishes with increasing dimen- 
sion. Therefore, no significant differences between the ex- 
ponents 0g''^ and v^'^ are observed in higher dimensions 
too. In the case of the three dimensional model sufficient 
data collapses (see Fig. ||) are obtained for the values 
9> = 1.53 ± 0.06, = 1.03 ± 0.03 and 9< = 1.54 ± 0.06, 
v< = 1.05 ±0.03. 

It is known that that the upper critical dimension of 
the CLG is 15c = 4 |l^. Thus one expects that a data- 
collapse of the corresponding curves is obtained if one uses 
the mean-field values of the exponent, i.e., i^y = 1 and 
9g = 3/2. In the latter case we assume that the time series 
of the order parameter can be mapped to a generalized 
random walk. Using these values a good data-collapse is 
obtained which is plotted in Fig.^. 
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5 Local persistence distributions 



We consider the local persistence distribution Pi{t) of the 
CLG model for D = 2,3,4,5. Compared to the global 
persistence fewer update steps are needed to determine 
the local persistence distribution for a similar degree of 
accuracy because the time evolution of all particles {N = 
L^p) can be used to measure Pi{t). Similar to the above 
analysis one varies again the exponents 0\ and i^i until one 
gets a data-collapse of different curves which correspond 
to different values of dp. The resulting data-collapses for 
D = 2,3,4,5 are plotted in Fig.^. The local persistence 
exponent decreases from di — 1.41 ± 0.02 for £> = 2 to 
6*1 = 1.19 ±0.03 for = 3. For > 4 good data-collapses 
are obtained for 9\ = 1, i.e., our results suggest that the 
mean- field value is (?i = 1. 

The values of the local correlation exponent fi de- 
creases from z/i = 1.10 ± 0.02 {D = 2) to vi = 1.04 ± 0.03 
{D — 3). In the case of D = 4 we used i^i = 1 to ob- 
tain the data-collapse. Surprisingly a significantly differ- 
ent value fi = 1.08 ± 0.03 has to be used to produce a 
data-collapse for = 5 (see Fig.^). But in the case of the 
five-dimensional system the simulations are limited by the 
available system sizes {L < 16). Therefore, it is possible 
that the above value describes a certain transient regime 
instead of the real critical scaling behavior. Therefore, we 
think that v\ — l describes the actual mean-field behavior. 

In the following we check the scaling relation Eq. (||) 
which connects the exponents /3, v\ and 9\ in the steady 
state. For D = 2 our numerically obtained result ia{2 — 
0i) = 0.649 agrees with the result /? = 0.637 ± 0.009 
from 0. In the case of L» = 3 we get iyi{2 - Oi) = 0.842 
which agrees again with (3 — 0.837 ± 0.015 obtained from 
the same reference. Finally the values vy = 1, 9i = 1, and 
(3—1 |lj,|l3 fulfill the scaling relation exactly which con- 
firms the above conclusion that i^i = 1 is valid for D > A. 

It is worth considering the average number of activa- 
tions within a certain time (riact ~ in detail. Using 
the Eqs. one sees that n^ct decreases with the dimen- 
sion, i.e., nf^l > for Sp^O because /3<°> < 
This suggests that the probability to re-activate a given 
site decreases with increasing dimension. This result could 
be surprising at first sight since one would expect that due 
to the increasing number of next neighbors in higher di- 
mensions it is more likely to re-activate a given lattice 
site. But the increasing number of next neighbors results 
in a reduced re-activation probability. We guess that this 
behavior is caused by a reduced return probability of ac- 
tivations, similar to a random walk. The probability that 
a symmetric random walk ever returns to its origin is cer- 
tain for D = 1 and D — 2 but uncertain for I? = 3 (see 
for instance (l8|). In the case of the CLG model one can 
assume that close to criticality only few activations diffuse 
through the system. Increasing the dimension it becomes 
less likely that these activations return to their origins. 
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Fig. 3. Scaling analysis of the global persistence distribution 
{t) for various dimensions. With increasing dimensions the 
slope of the curves (i.e. the exponent ^g) decreases slightly. The 
curves of the three- and four-dimensional system are shifted in 
the lower-left direction. The corresponding data-collapses for 
Pg (t) are of similar accuracy (not shown) . 




t (P-Po)"' 

Fig. 4. Scaling analysis of the local persistence distribution 
Pi{t) for various dimensions. With increasing dimensions the 
slope of the curves (i.e. the exponent 9i) decreases. For D — 
3, 4, 5 the corresponding curves are shifted in the lower-left 
direction. 

6 Discussion 

Considering the results of the local and global persistence 
exponents (see Table||) one concludes that the CLG model 
is characterized by the reversed persistence behavior 9g > 
Oi. The usual behavior 6g < 9\ was observed in reversible 
equilibrium phase transitions and means that the global 
order parameter persists longer in a given phase than 
the corresponding local quantities. For instance it is more 
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likely in a magnetic system that a single spin flips than 
the magnetization. 

Our result and the similar one of suggest that crit- 
ical systems with absorbing states behaves generally dif- 
ferent from reversible equilibrium systems. Here, the or- 
der parameter fluctuates on smaller scales than the local 
sites. Approaching the transition point (dp —>■ 0) only a 
few sites are activated. Considering snapshots of the CLG 
model we observed that these activated sites are concen- 
trated around small regions of activation which " diffuses" 
very slowly through the system. Thus non-active particles 
persist for a long time in their state until they are re- 
activated by other active particles. But the order parame- 
ter fluctuates around its mean value within this time. This 
behavior is a common feature of absorbing phase transi- 
tions and we therefore think that systems exhibiting an 
absorbing phase transitions are generally characterized by 
the reversed persistence behavior 9g> 9\. 



Table 1. The critical exponents of the CLG model for various 
dimensions D. The symbol * denotes logarithmic corrections to 
the power-law behavior. The values for pc, /3, 7' are obtained 
from The values of are obtained from the scaling 

relation 7 = v^D-213 



D 


Pc 


13 


7 


i'± 








Oi 


2 


0.34494 


0.637 


0.384 


0.83 


1.21 


1.10 


1.58 


1.41 


3 


0.2179 


0.837 


0.18 


0.62 


1.04 


1.04 


1.54 


1.19 


4 


0.1571 


1* 


0* 


1/2 


1 


1 


1.5 


1 


5 


0.1230 


1 





1/2 


1 


1.075 


1.5 


1 
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